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Recent results from the HI and ZEUS collaborations on hard QCD processes in 
ep interactions are reviewed. The topics cover jet shapes, the structure of real and 
virtual photons, and the production of J/^P mesons and open charm. 



1 Introduction 

Heavy flavours and jets in the hadronic final state carry much information 
about the dynamics of the unobservable quarks and gluons. They thus provide 
access to the underlying hard partonic processes of high energy ep interactions. 
The aim is to test the description of the production processes by perturbative 
QCD, and to extract information on the parton distributions of the initial state 
hadron. 

At HERA, electrons (or positrons) of 27.5 GeV collide with 820 GeV pro- 
tons, the center-of-mass energy is y/s — 300 GeV. The bulk of the scattering 
reactions occur in the photoproduction regime where the momentum transfer 
Q 2 w 0, i.e. the exchanged photon is almost real; data in the Deep Inelastic 
Scattering (DIS) regime correspond to large 7 virtualities Q 2 >> 1 GeV 2 /c 2 . 
In both regimes, however, the "fp center-of-mass energies W are typically as 
high as w 100 - 200 GeV and open up the phase space for hard QCD processes. 
The HERA experiments HI and ZEUS thus extend the range of fixed-target 
photoproduction experiments by an order of magnitude. 

This review begins with results from the study of jets in photoproduction 
and DIS, the focus is on the structure of both real and virtual photons. In the 
second part, the production of bound and open charm states will be discussed 
in terms of perturbative QCD models. 



2 Jets in Photoproduction 

In Leading Order QCD, a photon can either interact directly with a parton 
in the proton ("direct" photoproduction, 75 — * qg and 73 — > qq), or it can 
fluctuate into a quark- antiquark pair, which may form a vector meson state 

"Invited talk given at the XVIIth International Conference on Physics in Collision, Bristol, 
UK, June 25-27, 1997 



1 



ZEUS 1994 



O 
II 



- -gluons 

— quarks 



• ZEUS Data 
PYTHIA 
res+dir 

gluons 

quarks 




Ef > 1 4 GeV 



(a) 



-1 -0.5 0.5 1 1.5 



V 



ZEUS Preliminary 



; : « ■ 



' i * 



* ZEUS DATA DIS NC 1 994 (e* p) 
■ OPAL DATA (e* e") 

A CDF DATA (pp) 
m DO DATA (p p) 



(b) 



Figure 1: Jet shape (fraction of transverse energy in a sub-cone of radius r) for a jet defined 
by a cone of radius R = 1, (a) in photoproduction, at fixed r = 0.5 as function of rapidity, 
(b) in DIS at HERA, as function of r, compared with LEP and Tevatron data. 



or may develop a quark-gluon structure without forming a hadronic bound 
state ( "anomalous " characteristics of the photon) . A parton from that state 
interacts with a parton from the proton ("resolved" photoproduction; most 
important are qq' — > qq' , qg — > qg and gg — > gg). Since the parton carries 
only a fraction x 7 of the photon's momentum, the "resolved" events at HERA 
exhibit a stronger boost into the proton ( "forward" ) direction, towards positive 
pseudo-rapidities ?/.]] 

2.1 Jet Shapes 

A necessary prerequisite for the interpretation of jet photoproduction data in 
partonic terms is being able to model the internal structure of hadron jets. Let 
the jet be defined by a cone algorithm which combines hadrons within a cone of 
radius R, measured in the (77, <fi) plane. A suitable quantity to describe the jet 
shape is ^(r), defined as the fraction of the jet's transverse energy contained 
in a sub-cone of radius r. *S?(r) increases towards 1 with r approaching R; the 
more collimated the jet, the steeper the rise and the higher ^ at fixed r. 

In Fig. |l|a ^(r = 0.5) as measured by ZEUS 111 is displayed as a function 
of the jet pseudo-rapidity ?y, for jets with transverse energy E± > 14 GeV 

6 The pseudo-rapidity is defined as r\ = — lntan(f5/2), where 9 is the polar angle measured 
with respect to the 2 axis, taken to be along the direction of the incident proton beam. 
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and R = 1. The result clearly shows that the jets are broader in^he forward 
region. Predictions of the QCD Monte Carlo program PYTHIAn are shown 
for comparison. The thick curve which gives the best description of the data, 
has been obtained with the inclusion of multi-parton interactions in "resolved" 
processes. The predictions are also shown separately for jets initiated by either 
quarks or gluons. Gluon jets are broader since they have higher color charge 
so that QCD radiation is more likely to occur. In the generator program, the 
average broadening with rj can be traced to an increased fraction of gluon- 
initiated jets, which in turn is due to the increased contribution of "resolved" 
events in the forward direction. 

It is interesting to note .that in DIS, at high Q 2 > 100 GeV 2 /c 2 , the jet 
shapes measured at HERA El agree very well with those observed at LEP 
(Fig. |l|b). This is expected in the quark parton model, where a quark is 
"kicked" out of the proton and produces the jet. The function W(r) exhibits a 
distinctly steeper rise than that observed in pp interactionsu, where predomi- 
nantly gluon initiated jets are produced. 

Jet shapes in photoproduction have also been calculated analytically using 
Next to Leading Order (NLO) QCD matrix elements In order to emulate 
effects of the experimental jet finding procedures, a parameter R sep has to 
be introduced into the theoretical calculation; R sep specifies the maximum 
distance of two partons to be merged into one jet. Given the freedom in this 
jet finding parameter, the measured jet shapes can be reproduced. 

2.2 Dijet angular distribution 

In the Leading Order (LO) picture, the kinematics of the partonic 2^2 
process can be fully reconstructed in dijet events from the jets' energies and 
angles. Such an analysis can reveal dynamical features of the strong interac- 
tions mediating the jet production. 

One defines an observable related to the momentum fraction of the parton 
in the photon, x® BS = J2j e ts E± ex p( — rf et ) l^yEe where E e is the e beam 
energy, and y = 1 — E 1 /E e can be obtained from the scattered electron or from 
all final state hadrons. For the scattering angle in the partonic center-of-mass 
system one has cos#* = tanh^ 6 * 2 — rf etl )/2. Although beyond LO only the 
sum of "direct" and "resolved" processes can be unambiguously defined, and 
the 2 — > 2 kinematics formulae do not apply exactly anymore, these quantities 
- being defined as hadronic observables - are meaningful at any order. 

The uncorrected x 7 distribution from ZEUS i is shown in Fig. ||a. It is 
qualitatively reproduced by the Monte Carlo calculations shown in the same 
figure. For direct processes, x° BS is expected to be close to 1. The data is thus 
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Figure 2: (a) Uncorrected x® BS distribution for dijet events, (b) Dijet angular distribution 
in the parton parton center-of-mass system. 



split into a "direct" sample (x 7 > 0.75) and a "resolved" sample (x 7 < 0.75). 

In QCD, the dominant Leading Order diagrams for direct processes involve 
quark (spin |) exchange, whereas for the "resolved" case, gluon exchange (spin 
1) dominates. The angular distributions are thus predicted to be different: 
~ |1 — cosO*^ 1 in the first case, and a steeper |1 — cosf?*|~ 2 behavior in 
the latter. This expectation is nicely confirmed by the ZEUS measurements 
shown in Fig ||b. The distributions reflect the spin of the exchange as well 
as the relative contribution.s-.of direct and resolved processes to each sample. 
Moreover, NLO calculationsQ exhibit the same pattern and also agree well with 
the data. 



2.3 Photon Structure 

The jet cross sections evidently depend on the parton distributions in the pho- 
ton. At HERA, the proton acts as hadronic probe, and the data are sensitive 
to both the quark and the gluon content of the photon. The measurements 
have provided the first determination oflhc gluon momentum distribution by 
HlB. The recent analysis presented herecl investigates the scale dependence of 
the parton distribution and confronts it to QCD predictions. 

HI has measured the double differential dijet cross section as a function of 
x® BS and Ej,, the second highest transverse jet energy squared, in the range 
0.1 < x~ ( < 1 and E\ > 8 GeV 2 . The data have been analyzed in terms of 
an effective parton distribution - i.e. not differentiating quarks and following 



4 



p T 2 [GeV 2 ] 

1 10 10 2 10 3 




Figure 3: (a) Scale dependence of the effective parton density /J, , of the photon, (b) 
Dependence of the photon structure function , measured in two photon reactions, on the 
scale given by the virtuality Q 2 of the probing photon. /J,, was divided by 5 [see text) and 
is shown as function of the scale given by the ppiiton transverse momenta p 2 E3- The lines 
represent calculations in the FKP model Ej with different parameter settings. 



the concept of Ref. til This is possible since the most important leading or- 
der matrix elements in "resolved" photoproduction have angular distributions 
similar in shape, they mainly differ in overall magnitude by the ratio of color 
factors. Consequently, the cross section can be described by means of a single 
effective subprocess with matrix element \M$es\ 2 ■ One defines effective parton 
distributions for photon and proton: 

feff(x,P 2 ±)= (<l( X 'Pi)+Q( X 'P±)) + 49{x,Pl) , (1) 

flavours 

where q and g denote the quark and gluon densities, respectively, x the momen- 
tum fraction and p\ the factorization scale of the process given by the parton's 
transverse momentum squared. The dijet cross section then factorizes to good 
approximation, symbolically: 

a ~ r eff (x 7 ,pl) ■ f eff (x p ,pl) ■ \M SES \ 2 . (2) 

The effective parton distribution f2^ extracted from the double differen- 
tial cross section is shown in Fig. ||a as function of the scale p\ for the range 
0.4 < x 1 < 0.7. The data clearly show a direct dependence on the scale; This 
positive scaling violation is remarkably different from the behavior of hadronic 
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structure functions, where the dependence on the scale enters only via the 
QCD evolution equations and, at large-momentum fractions, is negative. For 
comparison the GRV parameterizationlill of the photon structure is shown (full 
line). In a vector meson dominance picture the effective parton density would 
be proportional to a vector meson structure function; here the effective parton 
density of the pion is shown, scaled by the probability of the photon fluctuating 
into a vector meson, this is also shown for comparison as the dashed line. The 
behavior of the photon is a genuine QCD prediction. It originates from the 
inhomogeneous term in the Altarelli-Parisi evolution equations which is only 
present for the photon, due to its point-like coupling to quarks, and is referred 
to as its "anomalous component" . 

The HI result ties in nicely with the scaling behavior of the photon struc- 
ture function as determined in two-photon experiments, where a highly 
virtual photon probes the quark content of a real photon. QCD predicts the 
scaling to be F^ ~ \nQ 2 / Aq CD . A similar dependence can be expected for 
f2ff at large cc 7 , sines here the gluon contribution to it is small according to 
the HI measurement u. From the GRV parameterization of quark and gluon 
distributions, the ratio f^ff/Fz is expected to be about 5 and approximately 
constant. Fig. ^dEI shows the HI result, scaled hy.this factor, overlaid onto a 
compilation of K7 results from e + e~ experiments!^. The HERA data are com- 
parable in precision, and they extend the kinematic range towards a maximum 
scale of p\ = 1250GeV 2 /c 2 . 



3 Jets in Deep Inelastic Scattering 

In the previous section, it was discussed how the structure of the photon can be 
"resolved" by partons inside the proton. The "resolution" is given by the hard 
scale of the partonic sub-process, e.g. the transverse energy squared, p\. On 
the other hand, deep inelastic ep scattering (DIS) is usually described in the 
quark parton model, where a virtual, point-like photon probes the structure of 
the proton, the "resolution" being given by the negative 4-momentum transfer 
Q 2 . In the kinematic regime of DIS, where Q 2 is low, but there is another hard 
scale p\ > Q 2 present, it is an open question whether this is still an adequate 
description. HERA offers the possibility to probe this interesting transition 
region between DIS and the photoproduction limit. 
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Figure 4: Dijet rate in DIS as function of Q 2 and XBj, compared to ARIADNE (full 
histogram), LEPTO (dashed) and DISENT (light dashed) calculations. 



3.1 Dijet Rates 

2+1 jet (dijet) events]] in QCD are due to 0(a s ) processes, namely radiation ol 
a hard gluon off the scattered quark (QCD-Compton), or boson gluon fusion 
(BGF). At low Q 2 < 100 GeV 2 /c 2 , BGF is strongly dominant. 

The dijet rate R2 measured by HlE_3 is shown in Fig. |as a function of Q 2 
and the Bjorken scaling variable XBj- The jets were identified in the hadronic 
center-of-mass system by a cone algorithm with R — 1 and were required to 
have E?p > 5 GeV. Various predictions are compared to the data. The Color 
Dipole Model (CDM) as implemented in tha ARIADNE Monte Carlo pro- 
gram describes the data well. The LEPTOliZI Monte Carlo contains the exact 
LO (C(a s )) matrix elements and parton showers in leading log approximation 



c The proton remnant jet is counted separately as 
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(MEPS). It falls short in describing the measured rate, in particular at lower 
Q 2 . The same holds for a full NLO QCD calculation (DISENT E3) where the 
jet algorithm is used to find jets from the final state partons. This discrepancy 
is unlikely to be an effect of missing hadronization effects, since these were 
found to be below 20% in the two Monte Carlo models. 

The failure of the NLO calculation to reproduce the dijet rate may be an 
indication that higher order corrections are important in the kinematic region 
probed here. The CDM also predicts additional energy flow (hadronic activity) 
in the photon fragmentation region of the jp center-of-mass system, increasing 
as Q 2 decreases, whereas in the parton shower model, this quantity remains 
small and constant. ,-Ihe Q 2 dependence of the HI data is reproduced by 
the CDM prediction E£l. It is possible to introduce additional processes with 
increased activity in this region into the cony£ntional QCD models by allowing 
or a partonic structure of the virtual photonESEI, analogous to the description 
of jet production by real photons. This will be investigated further in the 
future. 

3.2 Virtual Photon Structure 

The parton densities in the virtual photon are expected in QCD to be sup- 
pressed with increasing virtuality Q 2 ; for Q 2 > p 2 , the point-like component 
dominates. Fpi: example, according to the phenomenological ansatz of Drees 
and Godbole EJ, following the analysis of Borzumati and SchulerlEJ, the Q 2 - 
dependent suppression with respect to the densities in the real photon is pa- 
rameterized as 

UA^PlQ 2 )-k^P 2 t )-L with L = HiP llfj%^f )] (3) 

The parameter uj controls the onset of the suppression. For gluons the sup- 
pression is stronger than for quarks (~ L 2 ), since they have to be radiated off 
a parton that is already off-shell. 

A study into this direction is the measurement of inclusive differential jet 
cross, sections at lowest virtualities, in the region < Q 2 < 50 GeV 2 /c 2 by 
HlE3. The jets are identified here by means of a k t clustering algorithm with 
a resolution scale parameter E cut = 3 GeV. The measurement is expressed as 
a 7*p cross section, formally defined by dividing the ep cross section by a flux 
factor given in the Weizsacker- Williams approximation, thereby correcting for 
the trivial Q 2 dependence due to the photon propagator. The cross section as 
a function of the virtuality is shown in Fig. for different regions of transverse 
jet energy. The data is compared to predictions obtained with the HERWIG 
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Figure 5: 7*p jet cross section as function of Q in ranges of jet transverse energy. 



Monte Carlo program El, which includes in addition to the "direct", point- 
like process (shown as hatched histogram) a "resolved" contribution. The 
7* parton densities are modeled with the Drees Godbole ansatz and u> — 
lGeV 2 /c 2 . Only the sum of the two contributions (full histogram) describes 
the data. The "resolved" component is found to be important in the kinematic 
region where the jet transverse energy E 2 > Q 2 . 
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4 Charm I: J/vff Photoproduction 

In the conventional picture of ep interactions, charm is produced almost exclu- 
sively via boson gluon fusion. It is therefore considered as the "classical" way 
to probe the gluon content of the proton. This is one motivation to study the 
probing process more closely, but charm production is an interesting testing 
ground of perturbative QCD in its own right. 

4-1 Elastic J/\& Photoproduction 

Elastic J /ty production offers a very clean experimental signature: the proton 
remains intact, the vector meson is detected in one of its leptonic decay chan- 
nels, and apart from the ee or fi/j, pair there is nothing else in the detector. 
The energy dependence of the cross section is a characteristic feature of the 
production dynamics. In Regge phenomenology, it is related to a Pomeron tra- 
jectory a(t) = 1 + A + a't, where t is the 4-momentum transfer to the protoncJ. 
The cross section behaves like 

a~(W*) 2X ; (4) 

A is a universal constant for soft hadronic processes and of the order 0.08. ..0.1. 
On the other hand, it has been argued that due to the high mass of the J/^P 
mason, perturbative QCD should be applicable. In the model of Ryskin et 
al.c3, the dominant diagram contains a gluon ladder, and the cross section is 
therefore proportional to the square of the gluon density in the proton: 

a ~ [xg{x)f with x = Mj /x5 ,/W 2 . (5) 

Now the gluon density at low x behaves like xg(x) <~ x~ x ; for the W depen- 
dence again the form of Eq. || results, but A has a different meaning. The rise 
of the gluon density rather gives A values in the region 0.2. ..0.3. Thus, in this 
model a much steeper rise of the cross section with energy is expected. 

At HERA, the Regge picture is beautifully confirmed by measurements of 
the total photoproduction cross section, and by detailed investigations of the 
elastic production of light vector mesons - see Fig. ^| and Ref. Ea for a recent 
overview. However, the rise of the elastic J /fy cross section is found to be in 
clear disagreemenL-with a soft Pomeron picture. This is now seen with the 
HERA data aloneB, for example, ZEUS find 

4A = 0.92 ±0.14 ±0.10 . (6) 

The W dependence can be well described by the (leading order) perturbative 
QCD calculation in the Ryskin et al. model, the agreement with experimental 
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Figure 6: Total and elastic vector meson photoproduction cross section as function of CMS 
energy W JP . The solid lines are parameterizations based on Regge phenomenology. 

data is particularly good when the MRS (A') El parameterization of the gluon 
density in the proton is used. The sensitivity to the gluon density is temptingly 
high in this model. However, recently there have also been suggestions to 
explain the_ateep rise of the J cross section within a modified soft Pomeron 
frameworkEj. 
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4-2 Inelastic J/ty Photoproduction 

Inelastic J production is described in the Color Singlet Model (CSM) 
as production of a cc pair via boson gluon fusion. Another (relatively hard) 
gluon is radiated off one of the heavy quarks, in order to transform the pair 
into the color-neutral singlet state forming the J meson. Experimentally, 
the process can be separated by means of the elasticity variable z that can 
be calculated from the longitudinal momenta of the J /ip products and of all 
final state particles: z = (E — p z )ji^/{E — p 2 ) a ll- It represents, in the proton 
rest frame, the fraction of the photon's energy transferred to the J/^S/ meson. 
Diffractive production leads to z values close to 1, "resolved" contributions can 
be suppressed by a lower cut on z. _ 

The differential cross section, measured by HI and ZEUScH as function 
of the J/~$> transversajnomentum squared is shown in Fig. ^a and compared 
to QCD calculations^ in the CSM, in LO and NLO. Next to Leading Order 
corrections are substantial and clearly needed to reproduce the data. The pa- 
rameters in the theoretical calculation (charm mass m c = 1.4 GeV/c 2 , renor- 
malization and factorization scale p — \/2m c , and Aqcd = 200 MeV) are 
slightly "stretched" ; other choices tend to give a lower cross section normal- 
ization and would leave room for additional production mechanisms, e.g. via 
fragmentation l3 or through color octet states. 

Recently, the role of color octet contributions has been intensively dis- 
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cussed, as a possible explanation for thejonexpectedly high J/ty and ty' pro- 



duction rates measur- 
non-relativistic QCD 



at the Tevatron el In the factorization approach of 
J/ty mesons can also be produced via color-octet cc 
configurations. The transition to the singlet final state is treated as a soft pro- 
cess described by non-perturbative matrix elements, which must be obtained 
from experimental data. They have been extracted£3 from the CDF J/\& rates 
and yield a prediction for the HERA regime that can be confronted to data. 
The cross section measured by HI and ZEUS is shown as a function of the 
elasticity z in Fig. Qb. The color octet contributions are expected to enhance 



the cross section in the high z regionEj. No indication for such an enhancement 
can be seen in the data yet. However, there are considerable theoretical uncer- 
tainties in the prediction. As an example, the band in the figure indicates the 
range of uncertainty due to an effective-transverse momentum of the incoming 
partons, due to initial state radiationE3. 



5 Charm II: D* Production 

Open charm at HERA has mainly been reconstructed in the fragmentation 
and decay chain c — * D* + — ► D°ir + — > (K ~n + )n + exploiting the clean exper- 
imental signal in the mass difference distribution of m{K~ n + ) — m(K~ir + ). 

5.1 D* Production by Real Photons 

The total open charm photoproduction cross section has been measured at 
HERAl3 at different photon proton center-of-mass energies up to W « 230 
GeV and found to be an order of magnitude higher than at fixed target energies 
(Fig. H|). Together with a new point at "low" W w 90 GeV the data start to 
provide some lever arm on the slope of the cross section versus energy, which is 
related to the x dependence of the gluon density in the proton. However, there 
are large factors for the extrapolation (in rj and mainly p t ) from the visible to 
the total cross section involved, which depend on the parton densities in the 
proton and the photon. It is thus not straightforward to interpret the total 
cross section data in these terms. 

Currently there exist two approaches to calculate differential distribu- 
tions for_D* photoproduction in NLO QCD. In the so-called "massive" charm 
scheme Ell, charm quarks are treated as massive particles which are only gen- 
erated perturbatively in the final state, whereas the 3 lightest flavours (and 
gluons) are the only active partons in the initial state proton and photon. This 
scheme is valid as long as pt is not too large with respect to the charm mass, 
and indispensable for the calculation of the total cross section. In contrast, 
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in the "massless" scheme EUo, charm is also an active flavour in the parton 
distributions of the proton and the photon. This scheme should work particu- 
larly well for large p t 3> m c ; here, scale-dependent perturbative fragmentation 
functions are needed to control divergences. The "massless" scheme predicts 
a large (> 50%) contribution of "resolved" processes of the type eg — > eg initi- 
ated by a c quark from the photon.^ For comparison, the "resolved" component 
due to due to gg — > cc, present in both approaches, is only of order 10%. 

Differential cross section^as function of (pseudo-)rapidity and pt have been 
measured by HI and ZEUSL3, with lower p t cuts at 2.5 or 3 GeV/c. The bulk 
of the data thus lies in a region where pt is neither particularly "small" nor 
"large". In Fig. ^a), the ZEUS results are compared to NLO QCD calcula- 
tions in the two schemes. Within the - still large - experimental errors, the 
shape of the distribution is reproduced in both approaches. The normaliza- 
tion prefers the "massless' calculation, which has been performed using the 
Peterson-Zerwas form of the fragmentation function Cj with a fragmentation 
parameter e c = 0.064 as it was obtained from fits to e + e~ data within the same 
framework. Within the "massive" scheme, such an analysis has not yet been 
performed; the parameter e c = 0.060 is taken fmm a leading order analysis Cj 
of low energy e + e~ data. It has been argued c3 that rather a smaller value 

d This has no equivalent in the "massive" scheme, but one should bear in mind that only 
the sum of "direct" and "resolved" contributions is unambiguously defined at NLO. 
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e c ~ 0.020 would be consistent with the NLO calculation and reconcile it with 
the data. 

In order to obtain further information on the production process, ZEUS has 
performed a dijet analysis similar to that described in Sec. |2.2| and measured the 
uncorrected x® BS distribution for events containing a D* meson tj. The jets 
were identified using a cone algorithm or - shown here - a kt cluster algorithm 
and were required to have transverse energy above 4 GeV. The result is shown 
in Fig. ||b and compared to the predictions of the HERWIG leading order 
Monte Carlo program. The data at lower x*? B values need a substantial 
"resolved" photon contribution, which in the Monte Carlo is dominated by 
processes initiated by charm quarks from the photon. 

5.2 D* Production by Virtual Photons 

In deep inelastic scattering, due to the recoil of the scattered electron, D* 
production can be measured in the full range down to p t = in the hadronic 
center-of-mass system, although a lower cut on the transverse momentum in the 
laboratory frame is usually necessary to suppress combinatorial background. 
The data is then dominated by pt ~ m Cl so that one can expect a calculation 
in the "massive" scheme with 3 active flavours to provide a good description 
of the process in this regime. Fig. |l^ demonstrates that this is indeed the 
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Figure 10: Differential cross sections for deep inelastic D* production. 



case: the NLO calculationJay Harris and Smith performed in this scheme, 
agrees well with the dataCJ in both magnitude and shape. One should also 
note that due to the smallness of the charm transverse momenta, p 2 <C Q 2 is 
in general fulfilled, so that the virtual photon can be considered as a point-like 
probe to a good approximation. Deep inelastic charm production is therefore 
a very promising process to extract information about the gluon content of the 
proton. 
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6 Conclusion 

Jets and heavy quarks carry the information from the partonic interactions 
in ep collisions. The new HERA results in this field as discussed here can be 
successfully described in the framework of perturbative QCD. Due to the large 
kinematic range at HERA, complementary concepts have to be applied. Charm 
quarks appear as massive products of boson gluon fusion - or as massless par- 
tons in the proton and even the photon. Virtual photons can be considered as 
point-like probes of the hadron structure, but they also reveal their nature as 
composite hadronic objects as it is well explored now for their real counterpart. 
It is this variety in the HERA phenomenology that holds the promise of pro- 
viding keys to new dimensions in our understanding of the strong interaction. 
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